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The cosmological model of the modified Chaplygin gas interacting with cold dark matter is studied. 
Our attention is focused on the final state of universe in the model. It turns out that there exists 
a stable scaling solution, which provides the possibility to alleviate the coincidence problem. In 
addition, we investigate the effect of the coupling constants ci and C2 on the dynamical evolution of 
this model from the statefinder viewpoint. It is found that the coupUng constants play a significant 
role during the dynamical evolution of the interacting MCG model. Furthermore, we can distinguish 
this interacting model from other dark energy models in the s — r plane. 



PACS numbers: 98.80.Cq, 98.80.-k 

I. INTRODUCTION 

There is increasing evidence that our universe is presentlyin a state of cosmic accelerating expansion |py] . This result 
has been confirmed with observations via Supernovae la [2[, Cosmic Microwave Background anisotropics Large 
Scale Structure formation [4], baryon oscillations fs'] and weak lensing etc. These observations strongly suggest 
that the universe is spatially flat and is dominated by an extra component with negative pressure, dubbed as dark 
energy In the recent years, various candidates of dark energy have been proposed. The simplest candidate is the 

cosmological constant [9] . But this scenario is plagued by the severe fine-tuning problem and the coincidence problem. 
Other possible forms of dark energy include quintessenc e [iq| . k-essence 11|, phantom p^], Born-Infeld scalars [l^. 



quintom [14], tachyon field [lH, holographic dark energy [ij^aiid so on. Additionally, the conjecture that dark energy 
and dark matter can be unified by using the so-called Chaplygin gas (CG) obeying an exotic equation of state (EoS) 

has been investigated in several literatures [TtI]. Recently, the CG model was generalized to some possible forms [Tsj . 
As an alternative model, the EoS of the generalized Chaplygin gas (GCG) [19] reads 

P9 = (2) 

where < a < 1. It is clear that the case a = 1 corresponds to the CG. Within the framework of Friedmann- 
Robertson- Walker (FRW) cosmology, this EoS leads, after inserted into the relativistic energy conservation equation, 
to an energy density evolving as 

where B is an integration constant. Hence, we see that at early time the energy density behaves as a dust-like matter, 
while at late time it behaves like a cosmological constant. So the GCG model can also be interpreted as an entangled 
mixture of dark matter and dark energy. This dual role is at the heart of the surprising properties of the GCG 
model. However, these so-called unified dark matter models have been ruled out because they produce oscillations or 
exponential blowup of the matter power spectrum inconsistent with observation [2^. But no observation so far rule 
out the possibility of the GCG as dark energy though it is disfavored as dark matter. In fact, the dynamical properties 
of the GCG have been studied in Ref. [21], which indicates that the EoS of the GCG may cross the so-called phantom 
divide w = —I'li there exists an interaction between the GCG and dark matter. Therefore, it is worthwhile to further 
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understand such unified dark matter models by studying another candidate for the generalization of the CG, referred 
to as the modified Chaplygin gas (MCG) [H, [i^, which is characterized by a simple EoS 

P9=M - (4) 

where A, B and a are constants and < a < 1. The attractive feature of this model is that the EoS looks like that 
of two fluids, one obeying a perfect EoS p = Ap and the other being the GCG. From Eq. (4), it is easy to see that 
the MCG reduces to the GCG if A = and to the perfect fluid if B = 0. Accordingly, the evolution of the energy 
density is given by 

where C is the constant of integration. Thus the MCG behaves as a radiation (when A = 1/3) or a dust-like matter 
(when A = 0) at early stage, while as a cosmological constant at later stage. 

As we all know, observations at the level of the solar system severely constrain non-gravitational interactions of 
baryons, namely, non- minimal coupling between dark energy and ordinary matter fluids is strongly restricted by the 



experimental tests in the solar systems [2^. However, since the nature of dark sectors remains unknown, it is possible 
to have non-gravitational interactions between dark energy and dark matter. In this paper, we study the dynamical 
evolution of the interacting MCG model by considering an interaction term between the MCG and cold dark matter. 
In our scenario, we find that there exists a stable scaling solution, which is characterized by a constant ratio of the 
energy densities of the MCG and cold dark matter. This provides the possibility to alleviate the coincidence problem. 
Moreover, we find that the final state is determined by the parameters of the MCG, a, A, and the coupling constants 
ci and C2. The latter two parameters represent the transfer strength between the MCG and dark matter. Interestingly, 
we find that the EoS of the MCG Wg tends to a constant, which is only determined by the coupling constants ci 
and C2. However, both the EoS of the total cosmic fluid w and the deceleration parameter q tend to —1, which are 
independent of the choice of values for the parameters. This indicates that the cosmic doomsday is avoided and the 
universe enters to a de Sitter phase and thus accelerates forever. Further, we investigate the effect of the coupling 
constants ci and C2 on the dynamical evolution of this model from the statefinder viewpoint. By our analysis, we see 
that the coupling cantants ci and C2 play a significant role during the dynamics of the interacting MCG model. It is 
also worthwhile to note that, we can distinguish this interacting model from other dark energy models in the s — r 
plane. 

In Sec. II, we study the dynamical evolution of the interacting MCG. Then we apply the statefinder diagnosis 
to the interacting modified Chaplygin gas model for various different parameters in Sec. III. The conclusions are 
summarized in Sec. IV. 



II. DYNAMICS ANALYSIS 



In our scenario, the universe contains the MCG pg (as dark energy), cold dark matter and baryonic matter p^j. 
For a spatially flat universe, the Friedmann equation is 

where H is the Hubble parameter, /^^ = SttG and p is the total energy density (nature units c = h = 1 is used 
throughout the paper). Then differentiating the above equation with respect to cosmic time t and using the total 
energy conservation equation 

p + 3i/(p + p)=0, (7) 
where p is the total pressure of the background fluid, we can get the Raychaudhuri equation 

H = -^{pg^Pg^ Pm^Ph)^ (8) 

in which pg represents the pressure of the MCG. 
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We postulate that the two dark sectors interact through the interaction term Q and the baryonic matter only 
interacts gravitationally with the dark sectors. Then the continuity equation is written as 

Ph + 3i/(p5) = 0, (9) 



Pm + ^H{p^) = Q, (10) 

Pg^3H{pg^Pg) = -Q. (11) 

Clearly, Q is the rate of the energy density exchange in the dark sectors and the sign of Q determines the direction of 
energy transfer. A positive Q corresponds to the transfer of energy from dark energy to dark matter, while a negative 
Q represents the other way round. Due to the unknown nature of dark sectors, there is as yet no basis in fundamental 
theory for a special coupling between two dark sectors. So the interaction term Q discussed currently has to be chosen 
in a phenomenological way [25[. One possible choice for the interaction term is Q = 3H{ciPm -\- C2Pg) [l^, where ci 
and Co are coupling constants. This form was first proposed in [l^ and it is a more general form than those found 
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25|, |28|, which can be obtained when ci = C2 = c, ci = or C2 = 0. 



A. Stability Analysis 

To analyze the evolution of the dynamical system, we introduce the following dimensionless variables: 

2 2 2 

_ ^ Pg _ ^ Pg _ ^ Pm /-. r)\ 

""=3^' ^=3^' ^=aff^- ^ ' 

Accordingly, the density of the baryonic matter is determined by the Friedmann constraint (6) as 

|^ = l-x-2;, (13) 

which implies that < x -\- z < 1 and < x, z < 1. Furthermore, using these variables, the EoSs of the MCG and the 
total cosmic fluid are respectively given by 



Pg_^y 

Pg x' 



(14) 



w = ^= , , =y. (15) 



P Pg ^ Pm^ Pb 

The sound velocity and the deceleration parameter respectively read 

cl = ^ = -a^ + {l + a)A, (16) 

Opg X 

'/ = -| = -l + ^(l + y)- (17) 

Note that the condition for acceleration is w < — | and the physically meaningful range of the sound velocity is 
< < 1. Using Eqs. (6)- (11), we can obtain the following autonomous system: 

X = -3[(1 + C2)x + ^ + ci2;] +3x(l +^), (18) 
y = -3[(l + C2)x + ^ + ciz][-a^ + (l + a)A]+3^(l + ^), (19) 

X 

z = -3[-C2X + (l-ci)z]+3^(l + 7/), (20) 

where the prime denotes a derivative with respect to iV = Ina. We set the current scale factor by ao = 1. Then 
the current value of N reads A/'o = 0. Setting x = y = z = 0, we can obtain the critical points (x*,y*,z*) of the 
autonomous system as follows: 



point (a): + x,), 1 - ^^^{1 + x,)), 
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Point 


Existence 


Eigenvalues 


Stability 


It;* 


Acceleration 


(a) 


ci < 0, but ci 7^ — 1 

— (1 — yJ — Ci) < C2 < (J 
(V^+y^)^ < A < 1 


Ai = 3(A - Cl + C2)Xs, 
A2 = 2 ~ + j (1 + j , 

A3 - 3(1 + a) [1 + (^-^1+^2) (1 + ^^)] 


Unstable 


A — ci-\-C2 /I 1 ^ \ 
^ (1 + Xsj 


INo 


(&) 


ci < 0, but ci 7^ — 1 

-(1 - ^f^f < C2 < 

(V^+y^)^ <A<1 


Ai — -3(A - Cl + C2)Xs, 
A2 = |(A-Cl+C2)(l-Xs), 

A3 - 3(1 + a)[l + I^^^±£a) (1 _ ^^)] 


Saddle 




No 


(c) 


ci < 1, C2 > or 

Cl > 1,C2 < 

q; 1 — C1+C2 ^ 
1 + a 1 — Cl — 


Ai = -3, 

A2,3 = 2(cf-l){ C1+C2) 

-(1 - ci)[2 - Cl + C2 + (1 + a)A] ± \s} 


Stable if 
Cl < 1, C2 > and 

Saddle if 
Cl > 1 , C2 < and 

_^l-ci+C2 < ^ 
l + a 1 — Cl — 

l + a ^ 1 — ct ^ 


-1 


All 

a, A, Cl and C2 



TABLE I: The properties of the critical points for the interacting modified Chaplygin gas model. Here, the parameter As is 
defined in Eq. (25). 



• point (6): (^^1^(1 -X,), ^^-i±-^{l-x,), 1 - - x,)), 



• point (c): (,^, -1, j^). 
Here the parameter Xg is defined by 



4Aci 



{A-Ci+C2)^' 



(21) 



Thus, we can constrain the parameters in the model under the physically meaningful conditions, namely, < x* < 1 
and < c^^ < 1. Concretely, we analyze the existence conditions for these three points respectively: 

• points (a) and (b): Considering the physically meaningful range of = 1^^(1 ± Xg), we can obtain that 
< ^~^^+^^ (l±x,) < 1,^7^0, and moreover, Xs is real. Together with the constraint of the sound velocity c^^ = A, 
namely, < A < 1, we can conclude that the existence conditions of points (a) and {b) are ci < 0, but ci 7^ —1, 
-(1 - v"^)^ < C2 < and (v'^+ v"^)^ < A < 1. 

• point (c): According to the meaningful range of = ^rr^^q^? we know that < -[^^"^^^ < 1, and therefore, 
cl < 1, C2 > or cl > 1, C2 < 0. Furthermore, the constraint of the sound velocity c^^ implies that the range of A is 
~TTa ^~i-'ci^ ^ ^ < 1^(1 ~ ce "^~^J^^^ ). Thus, the existence conditions of point (c) are obtained. Results from our 
above analysis are concluded in Table I. 

To study the stability of the critical points for the autonomous system, we substitute linear perturbations x 
x^ + 5x^ y ^ -\- Sy and z ^ -\- Sz about the critical points into the autonomous system Eqs. (18)-(20). To 
first-order in the perturbations, one gets the following evolution equations of the linear perturbations: 

5x = 3{y^ — C2)Sx -\- 3{x^ — l)Sy — 3ci6z^ (22) 
6y = -3[a^(^+ci— ) + (l + C2)(l + a)A]fe 

+3[1 - (1 + a) A + (1 + C2)a + 2y, + a{2^ + ci — )]Sy - 3ci[-a— + (1 + a)A]Sz, (23) 

8z = 3c2Sx ^3z^8y ^3{y^ ^ ci)6z. (24) 

The three eigenvalues of the coefficient matrix of Eqs. (22)- (24) determine the stability of the critical points. We list 
the three eigenvalues for each point in Table I. Moreover, the condition for stability and acceleration (i.e., < — |) 
are also summarized in the same table. For convenience, we introduce the parameter 



A, = {[a(l-ci + C2) + (l+a)(l-ci)A]2+2(l-ci)(ci+C2)[a(l-ci+C2) + (l+a)(l-ci)A] + [(l-ci)(ci-C2)]'}i (25) 
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a 

FIG. 1: The stable regions in the {a, A) parameter space. The sohd and dotted hnes respectively denote the region for 
ci = C2 = (the modified Chaplygin gas model without interaction) or ci / 0,C2 = (coupUng between dark sectors only 
proportional to the energy density of dark matter) and ci = C2 = 0.2. In the region I (the region between the two same lines), 
point (c) is a stable accelerated attractor. The region II represents the region of solution without physical meaning. 




FIG. 2: The stable regions in the parameter space for the fixed ci or C2. The left two correspond to ci = (coupling between 
dark sectors only proportional to the energy density of dark energy) and ci = 0.4, respectively. The right is for C2 = 0.02. The 
region between the two colorful planes represents that point (c) is a stable accelerated attractor. The other regions correspond 
to the regions of solution without physical meaning. 

From Table I, we can clearly see that points (a) and {b) are not stable if they exist, and point (c) is stable under 
conditions ci < 1,C2 > and —j^^-^r^ < A < — ). Furthermore, the stable attractor is a 

scaling solution since the energy density of the MCG remains proportional to that of cold dark matter, = 
when C2 7^ 0. Thus, point (c) is an accelerated scaling solution that probably alleviates the coincidence problem. 
Additionally, for the stable attractor point (c), the location of the point depends only on the coupling constant ci and 
C2. The EoS of the total cosmic fluid reads = —1, which indicates that point (c) is a accelerated attractor {w:^ ^ ~|) 
and implies that the universe will enter to a de sitter phase and accelerate forever, and therefore, there is no singularity 
in the finite future for any parameters. However, the EoS of the MCG is expressed by Wg^ = —^-^z^^ = — 1 — jz^? 
that is to say, point (c) has a phantom equation of state {wg^ < —1) for ci < 1,C2 > 0. According to Eq. (17), we 
can obtain that = — 1 + |(1 + ?/*) = —1 for the stable attractor point (c), which also implies that the universe will 
enter to a de sitter phase and accelerate forever. 

B. Numerical Results 

In what follows, we numerically study the dynamical results from the autonomous system (18)-(20) to clearly 
confirm the complicated stability condition for the stable accelerated attractor point (c). In Figs.l and 2, we depict 
the parameter space for point (c) to be stable. We suppose that ci = C2 = c or ci 7^ 0, C2 = in Fig.l. From the 
figure, we see that the parameter space is independent of the coupling constant ci when C2 = 0, namely, the stable 
region of the modified Chaplygin gas without interaction (ci = C2 = 0) is the same as that of the interacting MCG, 
in which the interaction term is solely proportional to the energy density of dark matter (ci 7^ 0, C2 = 0). Fig. 2 shows 
that the stable regions in the parameter space for the fixed ci or C2. In the figure, the left two plots is respectively 
for ci = and ci = 0.4, meanwhile, the range of C2 is specially chosen as < C2 < 2 to present the results more 
transparently though C2 > is permitted. The right plot is shown by choosing C2 = 0.02, and also the range of ci is 
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FIG. 3: The evolution of the EoSs of the modified Chaplygin gas Wg for different values of various parameters. The top two 
plots correspond to the cases for different values of the modified Chaplygin gas parameters a or A, when ci = 0.2 and C2 = 0.4. 
In addition, we respectively fix the parameter A = 0.1 and a = 0.5 in these two plots. The bottom two plots represent the case 
for different values of the coupling constants ci or C2, when a = 0.5 and A = 0.1. Moreover, in these two plots, the parameter 
C2 and ci are fixed to be C2 = 0.4 and ci = 0.2 respectively. 

selected as — 1 < ci < 1 for clarity although its allowed range is ci < 1. 

In Figs. 3-4, we plot the evolution of the EoSs of the MCG Wg and the total cosmic fluid w for different parameters 
in the stable region to illustrate the evolutional feature of our universe. Meanwhile, the initial condition in these 
figures is taken to be xq = Qgo = 0.7, yo = —0.7 (we assume Wgo = — 1) and zq = ^mo = 0.26 to compare with 
the observation data. As example, we choose several sets of values for the parameters a, A, ci and C2 for clarity. 
From Fig. 3, we see that in the final state the EoS of the MCG, Wg^ tends to a constant, which is only determined 
by the coupling constants ci and C2. Furthermore, the final EoS of the MCG, Wg^^ is always below —1 for different 
parameters in the stable region, i.e., the MCG has a phantom equation of state in the final state. In Fig. 4, we find that 
the final EoS of the total cosmic fluid tends to —1, which is independent of any parameters in the stable region. 
This indicates that the cosmic doomsday is avoided and the universe accelerates forever. Note that the evolution of 
w is exotic when a set of parameters is chosen as a = 0.5, A = —0.1, ci = 0.2 and C2 = 0.4. In this case, the value 
of w is smaller than —1 at the high redshift, i.e., the interacting modified Chaplygin gas model behave as phantom 
{w < —1) at the not very early universe. 

As the most significant parameter from the viewpoint of observations, the deceleration parameter q is also discussed. 
We exhibit the evolution of the deceleration parameter q for various parameters in the stable region in Fig. 5. The 
sets of values for all parameters in this model are selected as those in Fig. 3-4. From the figure, we find that in the 
final state the deceleration parameter q tends to —1 for different values of various parameters. This further indicates 
the universe enters to a de Sitter phase in the final state. Moreover, varying any one of the MCG parameters a, 
A and the coupling constants ci, C2, and fixing the rest, we find that the larger the variational parameter is, the 
smaller the transition redshift Zf is, namely, the later the transition from the deceleration phase to the acceleration 
phase is. In addition, it is worthwhile to note that when we fix the parameters a = 0.5, ci = 0.2 and C2 = 0.4, and 
select A = —0.1, there are two transition redshift: the first one denotes the transition from the acceleration phase to 
the deceleration phase; the second one represents the transition point from the deceleration phase to the acceleration 
phase. This provides the possibility that the evolution of our universe from acceleration to deceleration and then from 
deceleration to acceleration. 
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-c, =0 
c = 0.4 



FIG. 4: The evolution of the EoSs of the total cosmic fluid w for different values of various parameters. In these four plots, we 
respectively select the same parameters as those in Fig. 3. 




FIG. 5: The evolution of the deceleration parameter q for different values of various parameters. Also, the parameters in these 
four plots are respectively chosen as those in Fig. 3 for clarity. 



III. STATEFINDER DIAGNOSIS 



In this section, we study the dynamics of the interacting modified Chaplygin gas model from the statefinder 
viewpoint. The so-called statefinder parameter was first introduced by Sahni et al. [HI in order to discriminate 
among more and more cosmological models. It is constructed from the scalar a and its derivatives up to the third 
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FIG. 6: The evolution of the statefinder for the interacting MCG model in the s — r plane with different coupling constants ci 
and C2, which is illustrated in the left one and the right two panels, respectively. To be clarity, we extract a part of the middle 
panel to be the right one. In these three plots, we fix the parameters a = 0.5 and A = —0.1. The solid points show the current 
values of the statefinder parameters (scTq). 



order. The statefinder pair {r, s} is defined as 



s = — (26) 



aH^' 3(g- 1/2)* 

Since different cosmological models exhibit qualitatively different trajectories of evolution in the s — r plane, the 
statefinder parameter is a good tool to distinguish cosmological models. It has a remarkable property for the basic 
spatially flat ACDM model and the matter dominated universe SCDM, i.e., the statefinder pair {r, s} for ACDM model 
takes the constant value (0, 1), while SCDM model corresponds to the fixed point (1, 1). We can clearly identify the 
"distance" from a given cosmological model to ACDM model in the s — r plane, such as the quintessence, the phantom, 
the Chaplygin gas, the holographic dark energy models, the interacting dark energy models, and so forth, which have 
been shown in the literatures [30|. 

Generally, according to the reexpression of the deceleration parameter q 

9 = -l-f, (27) 

we can also rewrite the statefinder pair {r, s} in terms of the Hubble parameter H and its first and second derivatives 
with respect to the cosmic time t, H and as 

r = 1 + 3-^ + -^, (28) 



3(^ + 1)' 



(29) 



Now we apply the statefinder parameter to the interacting MCG model. By using Eqs. (6)- (11), the statefiner 
parameter can be concretely expressed as 

r = l + ^[(l + C2)x + 7/ + ciz][-a- + (l + a)A], (30) 

_ [(1 + C2)x + ^ + ciz] [-gg + (1 + a) A] 
s - - . (31) 

The evolution of statefinder for the interacting MCG model in the s — r plane is plotted in Fig. 6. The arrows in the 
figure denote the evolution directions of the statefinder trajectories. In this figure, we also use the initial condition: 
xq — ^}go — 0.7, yo = —0.7 (we assume Wgo = — 1) and zq = Qrno = 0.26. Since the evolution trajectory in the s — r 
plane will be interrupted when y = {i.e.^w = 0), which is not our favorite case, we carefully choose the parameter 
A = —0.1 based on the model parameters selected before to avoid the interruption. This choice can be seen from 
Fig. 4, in which we can avoid the possibility of w = during its evolution only if A > based on the model parameters 
we selected. So in Fig. 6, we take A = —0.1 and a = 0.5 as example to investigate the effect of the coupling constants 
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ci and C2, and furthermore, to discriminate between the interacting MCG model and other dark energy models. From 
Fig. 6, we see that the evolution trajectories not only cross the ACDM fixed point but also end at the fixed point 
for various values of ci and C2. However, the r{s) curve could not traverse the SCDM fixed point for any case. In 
addition, in the left plot, the r{s) curves are very different with various ci and the fixed C2, especially, the current 
values of (5,r) are different for various values of ci. Also, we can clearly see the tremendous difference between the 
current values (5o,ro) and the ACDM fixed point, and the larger the value of ci is, the greater the "distance" from 
the current value (5o,ro) to the fixed point is. On the other hand, fixing the coupling constant ci and varying C2 as 
illustrated in the right two plots of Fig. 6, we find that the r{s) curves are very different from each other with various 
C2, and the current values of (s, r) are distinct for different values of C2. Also, we can clearly see the difference between 
the current values (so,^o) and the ACDM fixed point, and the larger the value of C2 is, the greater the "distance" 
from (so, To) to the fixed point is. Therefore, we can conclude that the statefinder diagnosis can not only discriminate 
the interacting MCG model with different coupling constant but also distinguish the interacting MCG model from 
other dark energy models. 

IV. CONCLUSIONS 

In previous sections, we have studied some physical properties of the interacting MCG model. By considering 
an interaction term between the MCG and cold dark matter, we study the dynamical evolution of this model and 
pay our attention to the final state of the universe. By our analysis, there exists a stable scaling solution, which is 
characterized by a constant ratio of the energy densities of the MCG and dark matter. This provides the possibility 
to alleviate the coincidence problem. Furthermore, we see that the final state is determined by the parameters of the 
MCG a, A and the coupling constants ci and C2. 

Interestingly, we find that the EoS of the MCG, Wg, tends to a constant, which is only determined by the coupling 
constants ci and C2. But both the EoS of the total cosmic fluid w and the deceleration parameter q tend to —1, which 
are independent of the choice of values for the parameters. This indicates that the cosmic doomsday is avoided and 
the universe enters to a de Sitter phase and thus accelerates forever. The transition from the deceleration phase to 
the acceleration phase also depends on the MCG parameters a, A and the coupling constants ci, C2. Varying one 
of these four parameters and fixing the rest, we find that the the larger the variational parameter is, the smaller 
the transition redshift Zt is, namely, the later the transition from the deceleration phase to the acceleration phase is. 
Further, we study the effect of the coupling constants Ci and C2 on the dynamical evolution of this model from the 
statefinder viewpoint. We clearly see that the coupling constants ci and C2 play a significant role during the dynamics 
of the interacting MCG model. The evolution trajectories of the statefinder for the interacting MCG model not only 
cross the ACDM fixed point but also end at the fixed point for various values of ci and C2. However, the r{s) curve 
could not traverse the SCDM fixed point for any case. Moreover, fixing one of the coupling constants ci and C2, and 
varying the other one, we find that the larger the variational parameter is, the greater the "distance" from the current 
value (so, vq) to the fixed point is. Thus, the statefinder parameter can discriminate the interacting MCG model with 
different coupling constant. It is also worthwhile to note that, we can distinguish this interacting model from other 
dark energy models in the s — r plane. 
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